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Abstract 

A recently developed lattice field theory approach to the statistical mechanics of charged 
polymers in electrolyte solutions [S. Tsonchev, R. D. Coalson, and A. Duncan, Phys. Rev. E 
60, 4257, (1999)] is generalized to the case where ground-state dominance in the polymer's 
Green's function does not apply. The full mean-field equations for the system are derived and 
are shown to possess a unique solution. The approach is applied to the problem of a charged 
Gaussian polymer chain confined to move within the region defined by two fused spheres. The 
failure of the notion of ground-state dominance under certain conditions even in the limit of 
large number of monomers is demonstrated. 
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1 Introduction 

The problem of partitioning of a polymer chain between cavities of different size 
has been of interest to researchers for some time [p], 0, |3|. Investigations of this 
phenomenon are motivated by the practical importance of techniques for separation 
of macromolecules according to their size, such as gel electrophoresis, size exclu- 
sion chromatography, membrane separation, filtration, etc. ||. All these methods 
rely on the different mobility of macromolecules of different size as they migrate 
through a porous network of random obstacles. Such networks can be modeled as 
a complex system of interconnected cavities and channels available to the polymer 
chain. Hence, understanding how a polymer partitions itself between such cavi- 
ties can lead to more efficient separation methods. Ultimately, one hopes to utilize 
the dependence of polymer partitioning on molecular properties, such as polymer 
length, electrical charge and electrolyte composition, to selectively separate chains 
according to their molecular weight. 

Recently, experimental and theoretical investigations f| § have explored the 
so called "entropic trapping" of polymer chains in large spherical voids in a gel, 
and have opened the door to potential new methods of macromolecular separation. 
"Entropic trapping" describes the process of preferential localization of a polymer 
chain inside a void of larger size than the typical channels and voids in a gel, due 
to the larger conformational entropy experienced by the polymer in the large void 
as compared to the one in the narrow channels and smaller voids of the gel. The 
trapping of the polymer in the cavity is characterized by the partition coefficient K, 
which in the case of a gel is defined as the ratio of the polymer concentrations inside 
and outside the cavity. 

In a recent letter [[/] we investigated the role of excluded volume interactions 
between monomers on the partitioning of a polymer chain between two connected 
spheres. In particular, we determined the dependence of the partition coefficient 
K (in this case, the ratio of the average number of monomers in each of the two 
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spheres) on the total number of monomers in the chain. The results, which are 
in accord with recent experiments || and computer simulations ||, showed that 
current theoretical understanding of polymer partitioning is incomplete. It was also 
shown that for certain kinds of systems the notion of ground-state dominance of the 
polymer's Green's function can fail even in the limit of large number of monomers. 

In this paper we investigate the role of electrostatic interactions on the partition- 
ing of a charged polymer chain between two confining spherical cavities of different 
size. We have carried out calculations using the Lattice Field Theory (LFT) ap- 
proach || , based on a formalism which has a number of antecedents in the literature 
[9-14]. 

In Section 2 we generalize the LFT approach to the situation where ground-state 
dominance of the polymer's Green's function is not assumed. Then in Section 3 we 
discuss the general shape of the total free energy functional of the system at the 
mean-field level, and prove that it has a unique minimum, thus guaranteeing the 
existence of a unique solution of the mean-field equations. Numerical results from 
calculations using three-dimensional (3D) LFT are presented in Section 4, and in 
Section 5 we summarize our conclusions. 
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2 Lattice Field Theory of Charged Polymer Chains in Elec- 
trolyte Solution 

We first extend the formalism presented in || to the general case where ground-state 
dominance of the polymer's Green's function does not apply. In Ref. |§ we showed 
that the full partition function of a charged polymer in an electrolyte solution with 
short-range monomer repulsion interactions can be written as a functional integral: 

Z = j D x (r^(r>^/* A ^-t/^^ , (1) 

where (3=1 /kT is the inverse temperature, e is the dielectric constant of the solution, 
e is the electron charge, A is a measure of the strength of the excluded volume 
interaction, x and oj are auxiliary fields, c±=e l3 ^ ± / \\ with /z± and A± being the 
chemical potentials and the thermal deBroglie wavelengths for the ions, respectively, 
and 

f f ds£ 2 (s)—ipeB f dsy(x(s))— i\ f dsu>(x(s)) 

Z Schr (x^) = / Dx(s)e ^ Jo PPJ J (2) 



with M being the total number of monomers in the chain, p the charge per monomer, 
and a p the Kuhn length. 

As before ||, the functional integral ([[]) can be rerouted through a complex 
saddle-point at x=— iXc and uj=— iu c , where Xc and u c are purely real, reducing 
the computation of Zg chr (x, oj) at the saddle-point to a conventional 3D Schrodinger 
Hamiltonian problem, that is, the computation of matrix elements of e~ HT , with 
Euclidean time extent of the evolution T=M and 

a 2 - 

H = V 2 + \u c (f) + ^ex c (r) . (3) 
o 

As usual, the equations determining the saddle-point configuration fields Xc u c are 
obtained by setting the variational derivative of the exponent in the full functional 
integral ([]]) to zero. In the general case of a polymer chain with free ends the polymer 
part of the partition function can be written: 

-ME„, 



7 j Sc hr = dxidxf^2^ n (xi)^ n (xf)e' 



n 

where E n is the n-th energy eigenvalue, 

A n = J drV n {r) , 



and 



(4) 



(5) 



(6) 



is the negative of the polymer contribution to the free energy. Thus the negative of 
the total free energy at the saddle-point becomes 



5 7T 



2 + ^ + c+e 0eXc + c.e"^ 
2 



+ F poi (x c ,w c ) . (7) 



As before [pfl, we have 



<5u; c (r) 



^pe|*„(f)| 2 , 



(9) 



Considering the variation of the potential energy V, 5V=j3pe5x c , we have 



(m|5V|n) 



(3pe E 



*m(r) 



E — E 



dr"m m {r")m n {r")8 Xc {r") 



so that 
Hence, 

and 



^Xc(r) Pe m w„ E n -E m 



(10) 



-ME r 



2{3peA n V n (r) ]T 



-ME„ 



(3peMA 2 n ^ n (r)e 



-MEr, 



(12) 



5F po ; 5 Xc (^ (^n e M£n ) 



j3pe 



E w , n ^gg^ e- M£ " -ME„ (A^ n (f)) 2 



(13) 
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Considering the fact that 

EA n ^ n A m ^/ m —ME n _ \ ■* A n ty n A m ty m _jvfi?„ 

n,m,m^n n m m,n,n^m m n 



and 



A n y n A m y m , MEn _ e _ ME ^ En ^E m _ M e -ME n ; (15) 

we can write (13) as an unconstrained double sum over states: 

^^^tibr^ -■ < 16 > 

In a similar fashion we obtain: 

r j-i v A„^„A m ^ m f -ME„ _ p — MEr, 

A v- a 9. — WB. • l i7 J 



As shown below, the quantity 

V AnjVAajgm ( p -ME n _ p -ME, 
2^n,m E n -E m V 



p(f)= - 



42 p-ME n v y 

is the total monomer density. Thus, we can finally write down the equations deter- 
mining the saddle point configuration fields Xc, u c : 

= -^V 2 x c (f) + c + e^^-c_e- /3e ^-pp(f) = 0, (19) 
pebXcXr) Aire 

1 8F 

\J^Jf) = ^) - p(r) = . (20) 
Using Eq. fl2Tj|) the auxiliary field uj c {r) can be eliminated, leaving the following 
pair of coupled nonlinear equations to describe, at mean-field level, the equilibrium 
properties of a charged polymer interacting with ions in an electrolyte solution: 

C V 2 X c{r) = c + e^ ex ^ - c_e-^ ( ^ - pp{r) , (21) 



4ne 



a 2 - 

-fV 2 * n (f) = Ap(r)*„(r) + [3pe. X c(r)Vn(r) - E n * n (r) . (22) 
o 

Eqs. (^1|) and (|22|) are generalizations of equations (20) and (21) in ||, as they 
involve the total monomer density, p(r), given by fllTf) , instead of M^^f), which 



would be appropriate only in the limit of ground-state dominance. The equations 
presented here apply for polymer chains of arbitrary length. Inclusion of single- 
particle potentials, which can be used to enforce exclusion regions for either the ions 
or monomers is straightforward ||. It is important to note that the parameters c± 
are exponentials of the chemical potentials p± for positively and negatively charged 
ions, the numbers of these ions must be fixed by suitably adjusting c± to satisfy the 
relations 

dlog(Z) 



n± = c±- 



dc~\ 



c± J e ±PeXc df. 



(23) 



Next we show that p(r) is in fact the monomer density. Starting from the average 



of the local monomer density, p n (f) fL5f , we have 
Pn(r) = Uir-Rn, 



JdR dR M (R M 


e -(M-n)H 




e -nH 


Rv) 


JdR dR 


M (Rm 


e -MH 


R 


•> 



Then the total density is 



M 



Pif) = Y,Pn(r) 

n=0 



M 



dnp n (r) 



Jo* dnfdR dR M (R M 


e ~(M-n)H 




e -nH 


Ro) 


JdR dRM ( 


Rm 


e -Mi 


i 


Ro) 



Invoking spectral decomposition 



r-2 



-mH 



n 



3=0 



E^(^)^(n)e 

3=0 



-mEj 



we obtain 



(24) 



(25) 



(26) 



dR^dRM (Rm 



-MH 



R ( 



E dR dR M ^ 3 {Ro)^ 3 {RNi)e 

3=0 J 



-ME, 



\2-ME, 



3=0 



(27) 
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Letting 



I in) = JdR dR M (R M \e- {M ~ n)(l 



r)(r 



e Ho 



^^•(it! M )^(r)e-( M - w )^ 

j=0 



E E AjA^^ke-^e-^-^ 

j=0 k=0 



(2? 



then 



pM oo oo „ m 

J dnl{n) = E AjA k ^^ k e- ME ' j dne~ n{ - E ^ E ^ 



where 



j=0 k=0 " 

oo oo 

EE^^/( M ;^»^ 

j=0 k=0 



for Ej ± E k 



(29) 



Ek—Ej 
-ME. 



Thus, all in all 



p(r) 



Ego J2 k=0 A 3 A k m^ k f{M- Ej, E k ) 
ET=oA 2 3 e- M ^ 



(30) 



thereby recovering Eq. (IS). It is easy to see that in the case of ground-state 
dominance, that is, if in the spectral decomposition formula (gfj) 

M(E 1 -E )»1, (31) 

we recover the familiar relationship 

p(f) = Mv^(f) , (32) 

and equations (|2lD and ( ^2]) become identical with equations (20) and (21) in ||. 

In concluding this section, we should mention that the functional F(x c , u c ) i n 
Eq. (^) represents a "mixed" free energy, since the partition function Z in Eq. ([!]) 
is canonical with respect to the polymer chain and grand canonical with respect to 
the mobile ions. The advantage of working with F is that, as will be shown in the 
next section, it has a unique minimum (corresponding to the solution of the mean 



field Eqs. (|21| ) and (0)), and thus can be used to guide a numerical search for 
the mean electrostatic and monomer density fields. Once the mean fields have been 
computed, the defining relation lnZ=F(x c , & c ) can De used to obtain free energies 
of various types. For example, the Helmholtz free energy A (corresponding to fixed 
numbers of monomers and impurity ions) is given by 

(3 A = n + In c + + n_ In c_ — F(xc <*> c ) • (33) 
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3 General Form of the Free Energy Functional 

In our previous work || we showed that in the case of ground-state dominance the 
negative of the free energy, F, is in fact convex, thus guaranteeing a unique solution 
to the field equations. Of course, convexity, although a sufficient condition, is not 
necessary for the existence of a unique local minimum of F. We shall see shortly 
that a weaker, but nevertheless sufficient condition for the existence of a unique 
minimum of F is the convexity of the functional e F . This latter convexity will be 
established below for arbitrary boundary conditions (open or closed) for the polymer 
chain, implying a unique minimum of F — a property essential for the stability and 
reliability of the numerical algorithms we employ to solve the mean-field equations 
of the theory. 

The convexity of e F implies the existence of a unique local minimum of F. This 
connection is due to the fact that local minima of F are inherited by e F , so that 
the presence of more than one local minimum of F would necessarily violate the 
convexity of e F . The exponential of the negative of the mean- field free energy given 
in Eq. (0) can be expressed as: 

Dx(s)e J 'G(uj c , £(s))H(xc, x(s)) , (34) 

where 

g(uj c ,x(s)) = J^dr-\Jdsu c (S(s)) ^ ( 35 ) 

H(Xc,x(s)) = e"^ i 'xcA X cdr+c+ J e^dr+c- J e-^dr-pe^ J ds X c{x(s)) _ /gg\ 

(This can also be seen by inspection of Eqs. (1,2) with x, u fixed at the saddle-point 
to — iXci —iu C) respectively.) Introducing the line distribution j(r)=Jds8(f—x(s)), 
Eqs. ( |35| , |36D become 

Q(lu c ,x(s)) = et/^-A/^rW)^ (37) 
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From Eq. ([34]), it is apparent that the exponential of F is a sum (over polymer 
configurations) of positively weighted products of functionals of the fields u c and 
Xc separately. Accordingly, the convexity of e F will be guaranteed once we can 
demonstrate the convexity of the functionals Q (as a function of uu c ) and Ti, (as a 
function of Xc) fo r any arbitrary fixed polymer configuration x(s). 

The second functional derivative of Q is found after a short calculation to be 

G{T,r ' ~ Siv c (r)Su c (f') 

= {xd(r- r ') + X 2 [u c (f) - j(r)] [u; c (f ') - j{r')]}Q . (39) 

A necessary and sufficient condition for a function W(f,r') to be positive semidefi- 
nite is: 

J dfj dr'f(r)W(r, r')f{r') > (40) 

for any function f(f). The positivity of the kernel Kg in Eq. (|3~9"D is apparent and 
can be demonstrated by verifying Eq. (|T0|): 

df J df'f(f)KG(r,f')f(f') 
= XQ j dff 2 (f) + X 2 gy drf(r) \u) c {f) -j(r)]} > 0. (41) 
The exponent of TC can be written: 

_ 8~~ / XcAx c dr + c + J e el3xc dr + c_ J e~ el3xc dr - pef3 J Xc(r)j(r)df 



~ J (xAXc + ^Xc(r)i(0) df+ c + J e e ^df+ c_ J , 



~ e ^dr 



8tt 



+c + y e eliXc dr + c„ J e 



e0Xc df 



+c+ / e e/3xc df+c_ / e^ e/3xc df, (42) 



where 



, _ 1 iTipe . 
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Hence, 



5K 

$Xc(f) 

and therefore 
K H (f,f') 



.Pi 

Air 



A (xc + ]) (r) + c + e(3e ePx ^ - cepe^*^ 



n 



(44) 



5 2 H 



^i A fxc + j) (r) + c + e(3e ePx ^ - c_e(3e- e/3x ^ 

4/T v ' 



$i 

47T 



n 



4/T 



+e 2 /? 2 <5(f - f ') (c+e e/3 * cW + c_e" e/3xc(r?) ) ft . 



(45) 



Property ( f40|) can be easily verified for the first and third terms of Eq. (f45| ). For 
the second term one finds: 



J dffdf'f(r) [-AfS(f- f ')] f(r')H = J dff(r) [-A/(r)] ft 



n>o. 



(46) 



This establishes the positivity of the kernel Kh, which together with the positivity 
of Kq proves the convexity of e F , and consequently the existence of unique minimum 
of F, which itself guarantees uniqueness of the solution of the mean-field equations. 
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Three-dimensional lattice field theory results for a charged 
polymer chain confined to two connected spheres of dif- 
ferent radii 



We have used a recently developed LFT approach || in order to obtain results of 
the solution of the mean-field equations ([H]) and (|22|) , for polymer chains confined 
to 3D cavities. In this section we present results of the calculation of the equilibrium 
monomer distribution of a charged polymer chain constrained to move inside the 
volume of two connected spheres |L6| . 

After rescaling according to f(r)— *fleXc(?)i ^N(r)— ^of *jv(^) 5 an d multiplying 
Eq. ( |2TD by af (a; being the lattice spacing), we solve the discretized version of 
equations (|2T| ) and fl22|) on a 3D lattice: 

aJ2^ftmffh = l+e fn -7_e _/ ™ - pp a (47) 

m 

%J2^r-trn^N,ra = ^X"Pn* N,n + Pfn^ N,n ~ E N ^ N>H , (48) 



r> / / j iif.il u j. v « / 1 1 

bar a. 



where 



* - 5*5. < 50 > 

and the wavefunctions are dimensionless and normalized according to 

= i; ( 51 ) 

n 

thus, the density sums to the total number of monomers, M. 

Equations (f|7]) and ( |48"D are solved simultaneously using the following relaxation 
procedure ||. First, the Schrodinger Eq. ( f48l) is solved for fa=0 and ignoring the 



nonlinear (monomer repulsion) potential term. The resulting ^a^'s and correspond- 
ing energy levels (wavefunctions and energy eigenvalues of a particle confined to 
a "box" consisting of two fused spheres) are used to calculate p^, then the Poisson- 
Boltzmann Eq. (f4"T|) is solved at each lattice point using a simple line minimization 
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procedure |L7|. The process is repeated and the coefficients 7± are updated after 
a few iterations until a predetermined accuracy is achieved. Then the resulting 
is used in Eq. (pESD, which is solved using the Lanczos method (TTSJ for a new set 
of ^ N,a s to be used in calculating an updated version of the monomer density p^. 
This density is then inputted into Eq. (|47D and a new version of fa is computed. 
For numerical stability, the updated fa inputted into Eq. (|48|) is obtained by adding 
a small fraction of the new (just obtained from Eq. (f47l)) to the old one (ob- 
tained from the previous iteration). The same "slow charging" procedure is used for 



updating p^ in the nonlinear potential term of the Schrodinger equation fl4"8|). 
The procedure is applied to the following two systems: 

(1) A polymer in a cavity consisting of two spheres of radii R\ and R2, carved 
out and sharing one common point on a lattice, the distance between the centers of 
the two spheres being R1+R2. 

(2) Same as (1), except that the spheres are now embedded in each other by one 
more lattice spacing, that is, the distance between their centers is now Ri+R 2 —a[. 

We have used the following parameters in relative units: i?i=1.0, i? 2 =0.8, a p =0.2, 
A=0.001, and the two spheres are carved inside a cube of 44 lattice points on each 
side with a; =0.1. In absolute units a p =5A. 

We have computed the partition coefficient K= (Mi) / (M 2 ) as a function of the 
total number of monomers in the system, M= (Mi) + (M 2 ), for varying monomer 
charge p and varying number of ions in the system. The results for system (1) are 
plotted in Figures 1 and 2, with n being the number of negative impurity ions in 
the system, while the number of positive ions is adjusted so that electroneutrality is 
preserved. Similar to what was found in the case of neutral polymers [|7) , we see here 
that for small M, InK increases sublinearly with M. This is followed by a turnover 
region, after which ln.fr decreases with M, and for very large M approaches a limit 
bounded from below by the log of the ratio of the volumes of the two spheres. 

In Figure 1 we show results for varying monomer charge p, keeping the number 



14 



500.0 



1000.0 

M 



1500.0 



Figure 1: The ratio \nK vs. M for system (1) for varying monomer charge p and fixed number of 
negative impurity ions n=600, which corresponds to a concentration C«0.75M. 




Figure 2: The ratio hiK vs. M for system (1) for varying number of negative impurity ions 
n and fixed monomer charge p. The corresponding concentrations are C(n=6)«0.75xlO -2 M, 
C(n=60)«0.75xlO- 1 M, and C(n=600)«0.75M, respectively 
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of impurity ions fixed. It is apparent that smaller monomer charge favors larger 
InK, thus making polymer separation easier. In Figure 2 we vary the number of 
impurity ions, showing that a large number of impurity ions leads to screening of 
the monomer charges from each other and thus to a behavior resembling that of a 
neutral chain. We observe an interesting feature in the case of small number of ions: 
the InK curve goes through a turnover, then through a minimum and a maximum, 
before approaching its asymptotic value at large M. 

In Figures 3 and 4 we show the analogous results to those of Figures 1 and 2, 
respectively, for system (2). The same basic behavior is observed, with the partition 
coefficient being smaller for system (2) than for system (1) under identical conditions 
of polymer length, monomer charge and impurity ion concentration. This is because 
of the wider conduit in system (2), which diminishes the isolation of one sphere from 
the other, thus enabling the polymer to move more freely between them. 



T ' 1 ■ r 




' I i i i i i i 

0.0 500.0 1000.0 1500.0 

M 



Figure 3: The ratio \nK vs. M for system (2) for varying monomer charge p and fixed number of 
negative impurity ions n=600, which corresponds to a concentration C«0.75M. 
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t ■ 1 ■ 1 ■ 1 ■ 1 ■ r 




0.70 1 ' 1 ■ 1 ' 1 ' 1 ■ 1 ' 1 ' 1 

0.0 200.0 400.0 600.0 800.0 1000.0 1200.0 1400.0 

M 

Figure 4: The ratio hiK vs. M for system (2) for varying number of negative impurity ions 
n and fixed monomer charge p. The corresponding concentrations are C(n=6)«0.75xlO _2 M, 
C(n=60)w0.75xlO" 1 Af, and C(n=60Q)«0.75M, respectively. 

These results are consistent with the ideas presented in our previous work ||, 
namely, that higher impurity ion (electrolyte) concentrations would favor better 
polymer separation between cavities of different size. 

Note that, as we found previously for an uncharged polymer chain @, in the case 
of system (1) ground-state dominance of the polymer's Green's function does not 
set in for any M. As we increase the total number of monomers, the first two energy 
levels come closer together and couple. Thus, even in the large M limit we must 
retain both these states in order to calculate an accurate monomer density p(r). 
In Figure 5 we plot A=exp \—M(E\—Eq)\ vs. M as an indicator of ground-state 
dominance (A— >0 as M^oo) for p=— 0.3 and n=600 for both systems (1) and (2). 
It is clear that in the case of system (1) ground-state dominance does not occur, 
while for system (2) it occurs, as would be expected for a polymer moving in a single 
cavity. 
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1 .00 




Figure 5: The number A= exp [— M (Ei — Eq)] as a function of M for system (1) (minimally- 
fused system) and system (2) (large-aperture system) for monomer charge p=— 0.3 and number of 
negative impurity ions n=600. 
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5 Conclusions 



We have extended the lattice field theory approach for the statistical mechanics of 
charged polymers in electrolyte solutions || to the case where ground-state domi- 
nance fails in the polymer's Green's function. At the mean-field level all thermo- 
dynamic properties of the system are obtained from the solution of two coupled 
nonlinear equations. These equations involve the full monomer density, which in 
general contains contributions from excited states in the spectral decomposition of 
the polymer's Green's function. We have also discussed the general shape of the neg- 
ative of the system's total free energy functional, and have shown that it possesses 
a single minimum, thus guaranteeing a unique solution of the mean-field equations. 

We have used this approach to calculate the equilibrium partition coefficient K 
of a Gaussian polymer chain in a system of two spheres connected by a narrow 
aperture, and have observed essentially the same generic behavior seen in the case 
of excluded volume interactions only 0]: the log of the partition coefficient increases 
sublinearly with the number of monomers in the chain, M, for small M, then it goes 
through a turnover region with a maximum, after which it decreases to an asymptotic 
value bounded from below by the log of the ratio of the volumes of the two spheres. 
Increasing the monomer charge makes this behavior less pronounced and reduces the 
maximum in the lni^ vs. M curve, while increasing the impurity ion concentration 
leads to screening of the monomer charges, and makes the lni^ vs. M curve similar 
to the one in the case of a neutral polymer with excluded volume interactions. This 
supports our previous contention || that higher impurity concentrations would lead 
to better polymer separation between cavities of different size, such as the cavities 
in a gel used to observe the "entropic trapping" phenomenon || . 

The present work demonstrates the failure of the notion of ground-state domi- 
nance in the case of a very narrow conduit between the two spheres. For such a 
system, even in the large M limit the first excited state of the polymer's Green's 
function must be retained in order to obtain an accurate mean-field solution to the 
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problem. 

Of course, the accuracy of the venerable mean-field approximation |T5|, |H| for the 



systems considered here remains an outstanding issue. In the case of electrically neu- 
tral polymers, Monte-Carlo simulations may provide valuable benchmarks against 
which mean-field predictions can be compared. Such computations are currently 
under way. 
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